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Abstract 

Horizontally  traveling  waves  in  the  ionosphere,  occurring  naturally  and  from 
nuclear  detonations,  cause  nonvertical  reflections  and,  therefore,  abnormal  iono¬ 
spheric  recordings  at  stations  during  overhead  passage.  It  is  shown  how  an  electron 
density  cross-section  in  the  vertical  plane  through  an  ionospheric  wave  can  be  con¬ 
structed  from  a  single  station's  ionospheric  recordings,  provided  those  are  taken  at 
time  intervals  not  exceeding  5  minutes  and  provided  the  general  direction  of  travel 
is  known.  The  described  analysis  yields  also  an  approximate  value  for  the  velocity. 
Figures  6  and  7  show  examples  of  profiles  of  ionospheric  waves  from  two  different 
nuclear  detonations,  observed  at  great  distances.  The  ionospheric  wave  of  30  Octo¬ 
ber  1961  from  Novaya  Zemlya  appears  to  be  caused  by  a  gravity  wave,  as  may  be 
inferred  from  Table  3  which  gives  a  review  of  world-wide  observations. 
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Electron  Density  Profiles  of  Wavemotions  in 
the  Ionosphere  Caused  by  Nuclear  Detonations 


1.  INTRODUCTION 

Ionospheric  h'f  recordings,  especially  those  taken  at  arctic  stations,  frequently 
during  one  24-hour  period  show  echoes  that  cannot  be  explained  by  assuming  hori¬ 
zontally  stratified  layers.  It  has  been  obvious  that  some  of  the  echoes  are  not  strictly 
vertical  soundings  but  arise  from  large-scale  irregularities  near  the  zenith  of  the 
station.  From  comparing  consecutive  h'f  recordings  it  has  also  been  obvious  that 
many  of  those  natural  irregularities  travel  with  a  speed  of  the  order  of  sound  waves. 
The  determination  o<‘ dimensions,  speed,  and  direction  is  possible,  in  principle,  if 
one  utilizes  several  sounding  stations  whose  spacings  are  large  enough  to  obtain, 
for  triangulation,  measurable  differences  in  echo-ranges  at  each  station  and  whose 
spacings  are  small  enough  to  obtain  recordings  that  can  be  correlated. 

Such  stringent  experimental  requirements  for  investigating  natural  large-scale 
movements  in  the  ionosphere,  requiring  three  ionospheric  stations  at  close  distance 
(100  km),  are  considerably  eased  fer  the  case  of  man-made  disturbances.  In  the 
latter  case  the  direction  of  movement  at  a  particular  station  and  the  speed  of  travel 
may  be  deduced  approximately  from  the  geographical  location  and  the  time  where 
and  when  the  man-made  disturbance  was  introduced.  With  this  information  avail¬ 
able  the  ionograms  from  a  single  station  are  sufficient  to  obtain  a  vertical  cross- 
section  of  the  ionospheric  layers  along  the  direction  of  travel  at  a  time  when  the 
disturbance  passes  through.  How  this  can  be  done  is  outlined  below  and  is  demon¬ 
strated  by  two  examples.  One  is  the  wave  from  the  nuclear  detonation  at  Novaya 
Zemlya  on  30  October  1961,  as  observed  near  Norway;  the  other  example  utilizes 
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ionospheric  data  taken  at  Kusaie  Island,  South  Pacific,  after  a  detonation  above 
Bikini  Atoll. 

2.  THE  RAW  DATA 

The  data  on  30  October  1961  were  taken  from  a  USAF  KC-135  jet  aircraft  at 
67oN-10°E  while  on  a  heading  toward  magnetic  north  (Figure  1).  The  distance  to 
Novaya  Zemlya  was  approximately  1600  km.  The  yield  of  the  detonation  was  esti¬ 
mated  as  being  between  35  and  90  MT;  and  the  altitude  of  the  burst  was  estimated 
as  being  below  the  tropopause. 1  Burst  time  was  0833  GMT.  The  equipment  in  the 
aircraft  consisted  of  a  modified  Step  Frequency  Ionospheric  Recorder*  that  sweeps 
from  1  to  25  Me  in  100  kc  steps  within  a  4-minute  time  period.  Nominal  transmitter 
output  was  above  10  kw.  The  transmitting-receiving  antenna  was  of  the  traveling- 
wave  type;  three  wires  approximately  23  m  long  extended  from  a  common  ti epoint 
on  the  aircraft's  upper  fuselage  to  three  separate  points  on  the  vertical  stabilizer, 
each  wire  terminated  by  a  resistor.  The  antenna  pattern  favored  upward  radiation 
in  a  wide  cone  of  about  90°  for  most  frequencies.  Anter.na  efficiency  was  below  but 
approaching  that  of  the  conventional  delta  antenna  of  an  ionospheric  ground  station 
*  most  frequencies. 

The  sequence  of  ionograms  on  35-mm  film  used  in  this  analysis  is  reproduced 
in  Figures  2a  and  2b.  The  same  ionograms  retraced  and  supplemented,  eliminating 
extraordinary  component  and  multiples,  appear  in  Figure  3. 

Data  from  the  Bikini  detonation  were  taken  at  Kusaie  Island  740  km  south  of 
Bikini  (Figure  4). ^  The  equipment  was  an  Ionospheric  Recorder  Type  C3,  which 
swept  the  frequency  range  from  1  to  25  Me  within  30  sec,  using  a  standard  delta 
antenna.  Detonation  occurred  approximately  1  hour  before  local  sunrise.  The 
detonation  above  Bikini  Atoll  was  below  the  tropopause  and  equivalent  to  several 
megatons  of  TOT. 

Figures  5a  and  5b  show  the  sequence  of  ionog-ams  used  in  this  analysis.  This 
sequence  covers  the  first  disturbance  arriving  at  Kusaie.  (A  second  disturbance, 
not  covered  in  this  analysis,  arrived  overhead  at  Kusaie  at  H+56  minutes. ) 


*  Manufactured  by  Philips  Electronics  Industries  Ltd,  Toronto,  Canada 
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Sequence  of  lonograms  taken  on  Aircraft  on  30  October  1961 


--S2 


Figure  2a, 
0920  GMT 


0850  to 


5 


Figure  2b.  Sequence  of  Ionograms  taken  on  Aircraft  on  20  October  1961,  0850  to 
0920  GMT 
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Figure  4.  Outline  Map  -  Bikini  Detonation  (in  both  geographical  and 
ordinates) 
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Figure  3a.  Sequence  of  lonograms  taken  at  Kusaie  Island  after  the  Bikini  Detonation 
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3.  THE  ELECTRON  DENSITY  PROFILES 


The  reproduced  sequences  of  ionograms  snow  that  the  single  echo  traces  near 
the  stan  of  the  sequence  develop  into  two  or  sometimes  mere  branches  (marked  I, 

II,  and  III  or.  Figure  3}.  This  means  that  for  certain  frequencies  reflections  occur 
at  more  than  one  distance.  This  indicates  clearly  the  existence  of  separate  ray 
paths  for  each  reflection  at  one  particular  frequency,  since  in  one  single  ray  path 
the  closest  reflection  surface  would  obscure  all  others  at  a  farther  hstance.  (Non¬ 
obscuring  partial  reflections,  which  produce  echoes  at  virtual  ranges  lacking  the 
feature  of  critical  frequency  and  associated  range  increase  (cusp),  are  not  seen  on 
the  ionograms  and  are  not  considered  here.  )  Each  branch  in  the  ionograms  is  the 
result  of  a  particular  set  of  ray  paths  (from  the  station  to  the  reflection  level  and 
back).  Each  set  may  consist  of  a  bundle  of  single  rays.  For  simplification,  how¬ 
ever,  it  has  been  assumed  that  each  br  nch  is  created  in  the  same  way  as  are  verti¬ 
cal  ionograms,  i.e.,  the  ray  paths  fo"  all  frequencies  of  a  branch  are  identical 
except  for  different  penetration  heights  into  the  layers.  With  this  simplification, 
reducing  the  cone  angle  of  the  bundle  of  ray  paths  tc  zero,  each  branch  of  the  iono- 
grsm  traces  can  be  reduced  from  virtual  to  true  height.  This  procedure  is  standard 
for  vertical  ionograms  and  is  done  with  the  notion  that  the  obtained  true  'height'  is 

actually  the  distance  to  the  layer  at  an  unknown  elevation  angle.  For  the  reduction 

3 

to  true  height  the  simple  method  of  Kelso  was  used  since  it  allows  rapid  hand  cal¬ 
culations  and  yields  sufficient  accuracy  compatible  with  that  of  the  other  steps  in 
the  analysis. 

Tables  1  and  2  contain  results  of  the  true  height  reduction,  referring  to  the 
ionograms  of  Figures  3  and  5,  respectively. 


4,  ESTIMATE  OF  VELOCITY 

The  next  step  in  this  analysis  is  the  postulation  that  the  irregularity  travels 
horizontally  as  a  whole  and  with  a  particular  velocity.  The  better  this  is  fulfilled 
and  the  better  the  velocity  is  estimated,  the  faster  the  analysis  described  under 
Section  5  leads  to  results. 

The  effect  of  the  Novaya  Zemlya  detonation  on  the  ionospheric  recordings  at  the 
aircraft  shows  up  first  at  G9C0  GMT  This  yields  an  average  speed  of  1600  km/ 27 min 
or  987  m/second.  On  the  other  hand,  other  widely-differing  velocity  values  rang¬ 
ing  from  330  tc  910  m/sec  from  ground  observation  were  measured  (Table  3  page  17). 
From  these  data  it  appeared  obvious  that  the  ionospheric  disturbance  had  started  with 
higher  velocity,  slowing  down  at  medium  distances,  and  that  the  aircraft  happened 
to  be  flying  at  a  transition  distance.  The  analysis,  therefore,  had  been  started  with 
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an  estimated  velocity  of  6  w  m/set  which,  however,  yielded  inconsistent  results. 
Several  lower  velocities  down  to  420  m/sec  were  then  assumed,  and  it  was  found 
that  most  consistent  results  were  obtained  with  a  velocity  of  500  m/second. 

The  arrival  at  Kusaie  of  the  ionospheric  disturbance  from  the  Bikini  detonation 
yielded  an  average  speed  of  514  m/second.  In  the  analysis,  however,  consistent  re¬ 
sults  were  obtained  with  velocities  centering  at  420  m/second. 


5.  CONSTRUCTION  OF  VERTICAL  PROFILES  OF  IONOSPHERIC  WAVES 

Figure  6  presents  a  two-dimensional  coordinate  system  showing  vertical  dis¬ 
tances  above  ground  and  horizontal  distances  from  an  arbitrary  zero-point.  The 
abscissa  also  shows  the  time  when  ionospheric  records  were  taken.  The  x*elative 
position  of  the  time  marks  to  the  kilometer  marks  on  the  abscissa  is  arbitrary, 
however,  the  value  9(distance)/9(time)  equals  500  m/sec  representing  the  speed 
postulated  above.  The  zero-km  mark  on  the  abscissa  may  be  considered  as  being 
fixed  with  the  horizontally-moving  ionospheric  disturbance;  the  ionospheric  station 
consequently  moves  along  the  zero  height  mark  from  left  to  right.  The  true  height 
values  (s  of  Table  1)  were  now  plotted  as  segments  of  circles,  with  that  particular 
radius  above  that  point  on  the  ground  which  is  marked  with  the  particular  time  when 
the  record  was  taken.  After  all  values  of  Table  1  were  plotted  in  this  manner,  a  first 
generation  of  contour  lines  of  the  plasma  frequency  f^  was  drawn  by  connecting  all 
segments  of  circles  of  the  same  plasma  frequency  by  smooth  curves.  Several  in¬ 
consistencies  remained,  especially  in  the  upper  levels  of  f^,  because  in  this  illustra¬ 
tion  each  ray  path  from  the  ionospheric  station  to  the  reflection  point  is  represented 
by  a  straight  line,  namely,  the  radius  of  the  particular  segment  to  the  point  where 
it  touches  the  contour  line. 

A  second  generation  of  contour  curves  was  then  obtained  by  replacing  those 
straight  ray  paths  of  obviously  incorrect  shape  with  curved  ones.  Samples  of  curved 
ray  paths  were  obtained  by  a  graphical  construction  according  to  Figure  8,  which  is 
based  on  Snell's  Law  and  which  approximates  the  refraction  process  by  a  step-wise 
procedure.  The  demand  for  sample  ray  paths  can  be  kept  to  a  minimum  if  one  is 
experienced  in  selecting  ray  paths  (choosing  the  right  departure  angle)  whose  re¬ 
flection  points  are  both  consistent  with  others  and  which  approximately  fulfill  the 
specular  reflection  requirement  (180°  return).  The  drawing  of  the  contour  lines  is 
then  a  matter  of  interpolation. 

A  third  generation  of  contour  lines  may  now  be  obtained,  commencing  at  the 
lowest  levels  of  uncertainty  and  proceeding  to  higher  levels,  by  calculating  the  virtual 
height  s'  =  ~  along  each  ray  path.  If  this  value  differs  by  A  s'  from  the  ex¬ 

perimental  s'  ('virtual  height'),  the  entire  procedure  might  be  repeated  with  new  values 


igure  6.  Cross-Section  of  Ionospheric  Wave  from  Novaya  Zemlya  Passing  over  Aircraft.  Dotted  c 
nes  are  not  measured  and  represent  suggested  supplementation. 
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Figure  8.  Graphical  Construction  of  Curved  Ray  Paths 

of  s  individually  adjusted  according  to  A  s'.  This  refinement  was  used  here  on 
only  a  few  points  as  a  check. 

The  described  construction  assumed  a  velocity  of  500  ±  80  m /second.  In  trials 
with  velocities  beyond  these  limits  certain  inconsistencies  could  not  be  removed. 

The  ionospheric  wave  over  Kusaie  was  constructed  in  the  same  way.  The  tabu¬ 
lated  values  of  Table  2  yield  contour  lines  as  seen  in  J  igure  7,  in  which  a  velocity 
of  420  ±  80  m/sec  is  assumed.  Assuming  velocities  b.yond  the  indicated  limits  led 
to  inconsistencies  during  the  construction  of  the  contour  lines. 

The  above  described  construction  of  contour  lines  may  be  applied  to  any  sequence 
of  ionospheric  recordings.  As  a  general  rule,  however,  it  is  felt  the  time  interval 
between  recordings  should  not  exceed  5  minutes. 

6.  DISCUSSION 

The  traveling  ionospheric  disturbance  from  the  30  October  1961  detonation  has 
been  observed  at  several  locations.  Table  3  shows  a  review  of  the  results  with 
velocities  deduced  from  the  various  data  and  assumptions.  The  values  bear  varying 
validity  margins  caused  by  the  uncertainties  of  distance  and  travel  time.  The  deter¬ 
mination  of  the  latter  is  affected  by  the  period  of  data  taking  and  by  the  uncertainty 
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of  the  time  of  origin.  If,  for  example,  the  ionospheric  wave  was  originated  by  the 
rising  fireball,  a  delay  after  detonation  of  as  much  as  5  minutes  should  be  assumed, 
raising  some  of  the  velocities  listed  to  even  higher  values. 

The  F-layer  disturbance  appears  to  have  started  with  a  velocity  exceeding  1000 
m/sec,  slowing  down  after  a  few  thousand  kilometers  of  travel  to  velocities  between 
330  and  430  m/second.  Simultaneously,  a  microbarographic  wave  propagated  to  far 
distances  at  lower  velocities  near  300  m/sec  as  a  separate  phenomenon.  The  F- 
layer  disturbance  consisted  of  a  sudden  onset  followed  by  consecutive  oscillations 
that  lasted  for  several  hours.  The  wave  appears  to  have  excluded  certain  directions 
of  propagation;  this  might  explain  the  discrepancy  in  the  round-the-world  velocity 

values  that  are  based  on  great  circle  paths.  The  observations  of  Stoffregen  and  of 

1  * 

Gardiner,  namely,  the  increase  of  the  amplitude  of  the  wave  with  distance  and/or 
with  decreasing  magnetic  latitude,  appear  significant. 

The  observed  wavelength  of  »  100  km  indicates  that  compressional  forces  play 

a  minor  role  in  this  wave  since  only  acoustic  waves  of  frequencies  >  0.  03  cps  or  of 

14 

wavelength  <  100  km  can  propagate.  Obayashi  suggests  that  the  ionospheric  dis¬ 
turbance  was  caused  by  a  gravitational  wave.  ^  Such  a  wave  would  have  large  ampli¬ 
tudes  at  great  heights.  However,  at  F-layer  maximum  and  above,  the  coupling  of 
motion  of  neutral  and  charged  particles  is  appreciable  and  imposes  an  attenuation 
on  the  gravity  wave.  By  this  effect  Obayashi  explains  why  the  largest  amplitude  of 
the  wave  is  observed  in  the  lower  F  layer. 

It  is  obvious  that,  under  the  influence  of  either  a  gravity  wave  or  a  compres¬ 
sional  wave,  the  motion  of  the  F-layer  plasma  will  not  coincide  with  that  of  the 
neutral  particles,  the  former  creating  dynamo  currents  and/or  electric  fields.  The 
cross-sections  shown  in  Figures  6  and  7  must  therefore  be  regarded  as  mere  ex¬ 
amples  of  a  variety  of  possible  forms.  It  must  be  expected  that  the  shapes  of  the 
ionization  contour  lines  will  vary  greatly,  depending  on  the  orientation  of  th^  geo¬ 
magnetic  vector  relative  to  the  direction  of  propagation. 
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